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Abstract 


Atherosclerosis is one of the most common diseases in the world. Med- 
ication with metal stents plays an important role in treating this disease. 
There are many models for delivering drugs from stents to the arterial 
wall. This paper presents a model that describes drug delivery from the 
stent coating layers to the arterial wall tissue. This model complements 
the previous models by considering the mechanical properties of the arte- 


rial wall tissue, which changes due to atherosclerosis and improves results 
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for designing stents. The stability behavior of the model is analyzed, and 
a number of numerical results are provided with explanations. A compar- 
ison between numerical and experimental results, which examine a more 


accurate match between the in vivo and in vitro, is shown. 


AMS subject classifications (2020): 65MO06. 


Keywords: Stent coating; Viscoelastic; Mathematical model; Numerical 


simulation. 


1 Introduction 


Atherosclerosis is a condition in which plaques made up of cholesterol, fat, 
cellular waste, calcium, and fibrin (a type of insoluble coagulation protein) 
form inside the arteries and block them [21] (see Figure 1). Methods used to 
reduce the severity of atherosclerosis include bypass grafting (BG), coronary 


skin intervention, and less invasive balloon angioplasty [12]. 


Figure 1: Atherosclerosis. Source: Intermountain Medical Center 
https: //healthcare-in-europe.com/en/news/potential-breakthrough-in-determining- 


who-s-at-risk. html 


Recently, though, a new procedure called stent-based delivery has been 
introduced that claims to help open blocked vessels and prevent restenosis 
[3, 28]. The bare metal stent is the first type of stent that is used as a scaffold 
to prevent restenosis of blood vessels (see Figure 2). Another kind of stent 
that is used to open blocked vessels and prevent restenosis is a drug-eluting 


stent [20]. In drug-eluting stents, the drug is gradually removed from the 


Iran. J. Numer. Anal. Optim., Vol. 14, No. 2, 2024, pp 475-499 


ATT Drug release from a two-layer stent coating considering the viscoelastic ... 


stent, which makes the drug enter the vessels with sufficient concentration in 
a timely manner and a biologically active state. This helps keep the arteries 
smooth and open and allows blood to flow well. Clinical evidence shows 
that drug effectiveness depends on release rate, drug diffusion rate, drug 
distribution in tissue, stent design, and so on. Therefore, it is necessary to 
study the relationship between stent design and the release and retention of 


eluted drugs. 


Figure 2: Stent. Source: Image by brgfx on Freepik.com 
https: //bjbas.springeropen.com/articles/10.1186/s43088-023-00382-9/figures/1 


Recently, researchers have studied different types of drug delivery based 
on stents, including spreading the drug from a metal stent surface as a Dirich- 
let condition [15, 22], single-layer coating of the drug on the metal wire 
[11, 16], and also multi-layer coating [14, 24]. For example, Sarifuddin and 
Mandal [24] have studied the benefits of replacing a single-layer coating with 
a multi-layer coating along with transfer, which includes better control over 
drug delivery. Multi-layer coating with different drug forms or release kinet- 
ics can provide better control over the drug transfer and treatment process. 

Atherosclerosis can affect the thickest arterial layer in advanced stages, 
and this effect can cause a significant change in the mechanical properties 
of the arterial wall tissue [10]. All of the above studies have examined drug 
transfer from stents to arterial wall tissue, regardless of the viscoelastic prop- 
erties of the arterial wall tissue. The model we present is based on a model 
presented in [24] and has the potential to help stent designers reduce adverse 
patient outcomes. The present study aims to investigate more precisely the 
drug delivery through a two-layer stent coating to the arterial wall tissue tak- 


ing into account the mechanical properties of the arterial wall tissue. More- 
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over, the model presented in this paper better agrees with the experimental 
results than the previous models. 

The paper is organized as follows: Section 2 presents a brief clinical per- 
spective of the work. Section 3 describes the model, including the model 
geometry, a mathematical description of the equations governing drug deliv- 
ery from the stent coating to the arterial wall tissue, and initial, boundary, 
and interface boundary conditions. The stability behavior of the model is an- 
alyzed in Section 4. Numerical simulations for drug delivery from stent layers 
to the arterial wall tissue are given in Section 5. Finally, some discussions 


and conclusions are mentioned in Sections 6 and 7. 


2 Clinical perspective 


Atherosclerosis is one of the most serious and common cardiovascular diseases 
in the artery wall. This disease causes blockage of the lumen and reduces 
blood flow. Surgical techniques, such as coronary artery BG and percu- 
taneous coronary intervention, are techniques used to reduce the severity of 
advanced atherosclerotic stenosis in the coronary artery. Balloon angioplasty 
was considered a less invasive intervention to open the blocked lumen. Clin- 
ical observations show that balloon angioplasty is associated with restenosis. 
The advent of bare metal stents has revolutionized balloon angioplasty by 
reducing the risk of restenosis. However, after several decades of using these 
types of stents, internal restenosis, that is, re-narrowing of the lumen, was 
seen again, which was a major defect. Therefore, the drug-eluting stent 
was invented to overcome the restenosis within the stent after implantation. 
Atherosclerosis mainly affects the intima, but the periphery, which is the 
thickest arterial layer, can be affected in advanced stages. This causes a 
significant change in the mechanical properties of the arterial wall. Several 
authors have studied mathematical models for drug delivery in cardiovas- 
cular tissues. Most of these studies focus on drug diffusion while ignoring 
the mechanical properties of the arterial wall. In this paper, we study the 
effect of the mechanical properties of the arterial wall, which has not been 
considered in the models that are used to deliver drugs from two-layer stents 


to the arterial wall. 
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3 Mathematical model 


This section describes the mathematical model of the problem, which consists 
of four parts: geometry of the model (Section 3.1), description of the drug in 
the stent coating layers (Section 3.2), description of the drug in the arterial 
wall tissue (Section 3.3), and interface boundary conditions (Section 3.4). 
The governing equations describe the release of the drug from stent coating 
layers S; (inner layer) and S$ (outer layer) into the arterial wall tissue V. The 
parameters used in the model are demonstrated in Table 1 with the relevant 


descriptions and units. 


Table 1: The values used in the model for the parameters. 


Variatle/Paxamevar Description Value Reference 
Co eeee maximum velocity (cm/s) 14 [24, 29] 
pt tissue density (g/cm) 1.06 5 
Lt viscosity (g/cms) 0.0072 5 
Vpite filtration velocity (cm/s) 5.8 x 107° [25] 
Pr permeability (cm?) 2x10 5 
Di tissue diffusivity (cm/s) 3.65 x 1078 6 
Di coating ($1) diffusivity (cm?/s) 1x 107° 6 
Dz coating (Sz) diffusivity (cm?/s) 1x 107° 6 
Dy viscoelastic diffusion coefficient (g/(cmsPa)) | 5x 1071! | Ref. model 
Bu tissue-binding capacity (mM) 1.3 13 
Kp constant (mM) 0.136 13 
Da Damkohler number 2700 2t 
T9 lumen radius (cm) 0.15 2, 30] 
ry thickness of coating (51) (em) 0.0025 24 
r2 thickness of coating ($2) (em) 0.0025 24 
T thickness of wall (cm) 0.05 2, 30] 
7 relaxation time (s) 0.5 Ref. model 
Te, relaxation time (s) 0.1 Ref. model 
K Young modulus (MPa) 0.0001 Ref. model 
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3.1 Geometry of the model 


This section describes the model’s geometry using a special stent shape with 
a half-embedded metal piece and a two-layer coating. To simplify the calcula- 
tions, a two-dimensional figure has been used instead of a three-dimensional 
one (see Figure 3). Also, S; and S2 show the inner and outer layers of the 
stent coating, and V shows the arterial wall tissue. The boundary between 
lumen and arterial wall tissue is called T;. [,, is a boundary between arte- 
rial wall tissue and other arterial wall tissue. The boundary between arterial 
wall tissue and perivascular wall is called [',. Moreover, I's,,y is an inter- 
face boundary between arterial wall tissue and the inner layer of the stent 
coating. The interface boundary between the inner and outer layers of the 
stent coating is called T's, 5,, and finally, the interface boundary between the 


outer layer of the stent coating and the metal part of the stent is called T's. 


Ta 
\ I's, 
Si T's,,S2 Pw 
So Ts 
ry Py 
(a) (b) 


Figure 3: A schematic of the geometric shape of the problem 


3.2 Description of the drug in the stent coating layers 


In this section, we will study the release of the drug from the stent coating 
layers. The drug is first loaded solidly into the stent coating. As the stent 
coating comes into contact with the fluid and the fluid enters the layers 
through their boundaries, the layer’s pores open. The fluid contacts the 
drug, and the drug begins to release. The drug released along the stent 


coating layers is described by the following non-dimensionalized equations: 
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Om, 
Ot 


=-V.J,,, in Sx (0,T;], m=1,2, (1) 


where Cm,m = 1,2 are the concentrations of the drug in the stent coat- 


ing (Sm,m = 1,2), respectively, J, 


Cm 


= —pi-Vem is the flux of the drug 
concentration in the stent coating (S,,,m = 1,2), and Ty is the final time. 
The parameters Pe,,,m = 1,2 are the Peclet numbers for the coating layers 


defined by Pe, = ols ,m=1,2, where D,,,m = 1,2 are the diffusion coef- 
ro 


ficient of the drug in coating layers S,,,m = 1,2, respectively, and Ro = @ 


and Up = Cus in which rp and Umax denote the luminal radius and the 


maximum velocity at the inlet, respectively. 


At time t = 0, we assume that the amount of drug released after contact- 


ing the loaded drug with the fluid in the stent coating is as follows: 


Cm(0) = 1, m= 1,2. (2) 


We assume a non-flux boundary condition on the Tg, that J.,.7 = 0 because 
the drug is impermeable to the metal part of the stent, where 77 is the exterior 
unit normal. The drug enters the inner stent coating via boundary condi- 
tions from the outer stent coating and then enters the arterial wall tissue 
through the inner stent coating. Section 3.4 describes the interface boundary 


conditions. 


3.3 Description of the drug in the arterial wall tissue 


The equations used to describe the distribution of the drug in the wall are 


as follows: 


Convection is caused by the flow of interstitial fluid within the porous 
arterial wall, and the continuity equation can be written in dimensionless form 
in two-dimensional form. The Brinkman equation provides more accurate 


results because it combines the linear momentum and mass conservation of 
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fluids in large pores and channels with the Darcy equation for regions with 
unresolved pores [7]. Therefore, the Brinkman equation is used as follows: 
1 
x +(v.V)v =-VP+ Fe ual +(Vv)7)- 
V.v=0 in V x (0,T;], 


1 
Vv, 
ReKy (3) 


where Re = a and Ky, = ra in which P, is the Darcy permeability. 
Here, v is the velocity vector for the tissue domain, p; is the density, juz is 
the viscosity, and P is the pressure. 

To complete the equation (3), the following boundary and initial con- 
ditions are given: We assume the initial condition v(0) = 0 and boundary 


conditions on the boundaries T’;,T, and T,, (see Figure 3), as follows: 


v=0 on Ty, 
V = (vz, vy) = (0, Vee) on Ty, (4) 
P=0 onTlag, 
where Viz is the filtration velocity of plasma. 
A drug that does not bind to tissue receptors is known as a free drug, the 
concentration which is indicated by c;. Its equation in the arterial wall tissue 


is as follows: 


Ocr Ob; . 
= Vide V 0, T b 5 
= .- SE in Vx OT) (5) 
where J., = —( oe Vo 4 ie Vo — vc) is the free drug mass flux, and o 


is the response of arterial wall tissue to strain caused by drug molecules. 


RoVo 
Dy, ? 


D, is the diffusion coefficient of the drug in the arterial wall tissue. Also, 


Moreover, Pe; is the Peclet number for tissue defined by Pe; = where 


D, = £24, where D, is the viscoelastic diffusion coefficient. Finally, b; is 


D 
the concentration of a drug that attaches to tissue receptors in the arterial 


wall and is known as a bound drug, which is described by the following 


equation: a 
a = kac:(Bur — 4) — kab, in V x (0,T;|, (6) 
where kg = Koto and kg = te in which K, and Kg are the aggregation 


and disintegration rate constants, respectively, where K, = 22s and Kg = 


BuT? 
KyKp. 
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The second sentence of J., in equation (5), which is BVO, describes the 
mechanical behavior of the arterial wall tissue. In other words, it describes the 
arterial wall as a barrier that prevents the drug from entering the wall. Using 
the same symbols used in [8] for stress and viscoelastic diffusion coefficient, 
the model from [4] 


in Ve OT, (7) 


is similar to the equation used by [19]. Also, « is the Young modulus, and + 


and 7, are the relaxation time. 


The initial and boundary conditions are taken into account to complete 
the equations (5) and (7). As in the initial time, there is no free drug and 


bound drug in the arterial wall, so we have 


cr(0) = b,(0) = 0. (8) 


Due to the fact that the free drug does not cross the arterial walls (T,,UT a) 
as well as the boundary between the lumen and the tissue of the arterial wall 


(T,), there is a boundary condition J.,.7 = 0 on the boundaries [; UT, UF a. 


3.4 Interface boundary conditions 


Aside from the initial and boundary conditions, Systems (1), (3), (5), and 
(7) are completed by the interface boundary conditions between the first and 
second layers of the stent coating (Is,,s,) and the second layer of the stent 


coating and the arterial wall (I's, v) as follows: 


Je, = —Agcg on V's, \5,, 
Jo. = —JI.,.v on Vg, \5,, 
Je, = Aic, on V5, v, 

Jo, 1 = —Je,.v on I's, v, 


where A, and Ap are partition coefficients. 
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4 Energy estimates 


In this section, we obtain an upper bound for the energy estimate, which is 


defined as follows: 


E(e(t)) = lle(e)|? + fl | Ve(s)||2ds + i le(s)ll2ng, eure, vy %: 


where c(t) = (c1(t), c2(t), %(t)) and || - || represent the usual norm in L?(Q) 
(Q = S,US2UV), which is induced by the corresponding inner product (-,-) 
[1]. After solving (7), we get the following equation: 


t ‘ —t —t 
c= aCe — yf er §Yer(s)ds + Ze e4(0) _ re +e7a(0). (10) 
0 


st 


Replacing (6) in (5) to get 


Oc 1 1 

= F a Vas D. Vo vcr) kact (Bur _ b:) + kaby. (11) 
The weak formulation of the differential problem (1) and (11) is defined as 
follows: Find c(t) € H1(S1) x H1(S2) x H'(V) such that 22(t) € L?(S;) x 


L?(S2) x E7(V), e(0) = (1,1,0), and 


1 
(5, (4), 4) a ¥ (— 5, Vein Vu) ac (Agce, ui)rs, 59 ~ (Aici, ui)rs,.v 


i=1,2,t 


1 
—(A2c2, ta)rs,.s, — (q-Vo — vee), Vue) + (Aner, Ue)rs,,v 


—(kac:(Bur = bi), Ut) + (kab, Ut); (12) 


where (u1, U2, ut) € (H1(S1) x H+(S2) x H1(V)). 
Replacing (10) in (12) by considering o(0) and c;(0) as a constant to get 


Oc 1 
(5,4) = ¥ (— pe, Vow Vu) + (Agéa; ti )rs, 2, — (Arei; M4) rs, v 
i=1,2,t ¢ 


K To 


—(Azea, u2)rs, 55 + D,7! = r 


oc) cx(s ut(s))ds 
)| (Ver(s), Vur(s))d 


KTo 
7) (Ver, Vue) + (vey, Vuz) + (Aici, Ui )T's,.v 
—(kace( Bur — bt), ue) + (Rabe, ut). (13) 
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(ci (t), co(t), (t)) is a solution to the variational 


Theorem 1. If c(t) = 
problem (13), then 


1 
E(c(t)) < ——E(c(0))e™, (14) 
where ¢ = fos is a time independent positive constant, with dmin = 
min{1, 2 Per? es Ces Dee €3),+2A1,2A2} and 
1 2 9 Ag 
Omax => max{(57yllvllz-(1~) = 2k,Bmu = 2ka6||be|| (L)), 26) Ao, 22 
K, ie Ay 
1— +), 263 Ai, 5 
Dor tp) 24a pa 


where €1, €2,€3 # 0 are positive constants and (p,- — p°% — 6) > 0. 


Proof. In (13), replace wi, ug, and uz with c1, co and c, respectively, to get 


lle)? =-Br- DL [WVeill? + (Arce, er)rs,.s. ~ Arllerlits, v 


* §=1,2,t 


NW] re 
| 


t 
K Te Lifgu: 
~Aalleals, og + Beg ~ 2) fet IVeu(s)Fas 


+ (wee, Vet) + (Arer, Ct) rs, .v 


~(haee(Bus — bt), cz) + (Kade, cr). (15) 


Using the Cauchy—Schwarz inequality and since kg = K pk, assuming 6 > 1 


such that c, > fa 
1 a Pi 
aac (t)|? < + fAalleallz,. So ae 4e2 lleallE.,, So 
4=1,2,t 


sin v ~Aalleallés, s, 


ane ga KT 
pet) f ere P vers) Pas+ 5 Ver? 


1 
Fe3|[Verll® + Fallvllzcocr~ylleell? + BAsleilBe # 


Ai 
re llcellts, y ~ Ka Barlleell? + Kabllbell pz) lleell?. 


Let €2 > 0 such that (p= — $= — 3) > 0 and (zezll¥llZ0 (2) — —kaBu + 


ka ||bel| r--(,°)) > 0. Then 
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1 1 2 1 KTo 2 2 2 
a ailcl oe Weil? + (Be poe allVerl? + Arlleallts, 


1 
+Aalleallrs, sy S (Gallvllix (r=) — kaBa + had |[belln~ (cy) lleell? 
2 


Ag 
+eAallealits, o. + Goals, 5 


t 
tpg By fet Miva) Pas 
oT T 0 
2 2 Aj 2 
+e Ail[eillp. +r de llellPs, v- (16) 
3 


Multiplying (16) by two and taking integration to get 


1 
dmin 


using the Gronwall lemma [9] to obtain (14). 


E(c(t)) < 


E(c(0)) +4 | E(e(s))ds, 


The estimation we obtained in Theorem 1 allows us to determine the 


uniqueness of the solution (13) and its stability for a limited time. 


5 Numerical simulations 


COMSOL 5.3 Multiphysics is utilized for meshing and numerical results. 
The arterial wall and coating layers are meshed, as shown in Figure 4. The 
following modules in the COMSOL are used: 


e Brinkman equations in order to calculate the velocity and pressure of 


the arterial wall; 


¢ Coefficient form partial differential equation (c) in order to calculate 


the drug concentration in the coating layers and arterial wall tissue. 


To make sure that the numerical results do not depend on the mesh size, 


we consider the following meshes: 


e Mesh 1: 1642 elements; 
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Figure 4: Computational meshes in the domain 


Table 2: Different errors for free drug concentration in arterial wall. 


h EO [? Hi 
0.001 1.3x10-%  154x10-" 3.81 x 190-1 
0.01 2.2x107!2 2114x1071! 5.88 x 10710 
0.1 3.6x 10712 3.56x 1071! 7.43 x 10719 


e Mesh 2: 3132 elements; 
e Mesh 3: 6162 elements; 
e Mesh 4: 12340 elements. 


The adequacy of the mesh density has been investigated using convergence 
tests. By considering two different meshes, the solutions’ convergence and the 
mesh’s independence have been investigated. To illustrate this claim more 
clearly, Figure 5 and Table 2 show the comparison among different meshes for 
the average drug concentration in arterial wall tissue. Of note, the reference 
solution is defined by hmax = 1 x 107+. According to this figure, the changes 
among the meshes are almost unrecognizable, suggesting that the numerical 
results do not depend on the mesh size. 

Based on Figure 6, the velocity vectors affect both the distribution and 


retention of the drug. This result indicates the boundary condition for the 
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0.01; 
0.009; 
0.008; 
0.007; 
0.006; 
0.005; 
0.004; 
0.003; 
0.002; 
0.001; 


Average concentration c_t 


0 1 2 3 4 5 
Time 


Figure 5: Average drug concentration (c;) in the arterial wall. 


interface boundary between the lumen and the arterial wall tissue. The 
velocity vector is indicated by an arrow. These results are consistent with 
the results in [28]. 


x10° 


p 


WwW 


Figure 6: Dissemination of the velocity vector 


The filtration velocity used in [24], which is derived from experimental 
results performed on animals [25], is similar to the amounts obtained in [17, 


23, 30]. Figure 7 shows the low effect of filtration velocity on the average 
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drug concentration in arterial wall tissue for 1.8 x 10~°(em/s) and 5.8 x 
10-® (cm/s). 


0.0106 


0.01 
0.009 
0.008} 
0.007 
0.006 
0,005} 
0.004 
0.003} 
0.002, 
0.001 


Average concentration ct 
Average concentration c_t 


0.0106 
0 1 2 3 4 5 0.88 0.89 0.9 0.91 0.92 
Time Time 


(a) (b) zoom of (a) around t=0.9 


Figure 7: The effect of filtration velocity (V+ ¢) on the average drug concentration c; in 


the arterial wall tissue. 


The release of drug concentration from the stent coating layers into the 
arterial wall tissue during 300 seconds has been demonstrated in Figure 8. As 
shown in Figure 8, the drug concentration enters from the inner layer of the 
stent coating due to the interface boundary condition on I's,,y. Over time, 
the drug concentration enters from the inner layer of the stent coating and 


increases in the arterial wall tissue. 


The local concentration of the free and bound drug, which is indicated by 
a red point (see Figure 9 (f)), is illustrated in Figure 9. According to these 
figures, the concentration of free drug increases for a certain period of time 
and then diminishes (Figure 9 (a) and (c)), although the bound drug goes 
up to reach a steady state (Figure 9 (b) and (d)). By comparing the rate of 
drug release from the inner and outer layers of the stent, it can be seen that 
a faster release leads to a higher concentration of the free and bounded drug 
in the arterial wall. For further explanation, we need to point out that the 
rapid diffusion in the inner layer compared to the outer layer has a higher 
concentration in the arterial wall tissue, which is around 0.16 in the inner 
layer (Figure 9 (a)) and around 0.12 in the outer layer (Figure 9 (c)). This is 
because the inner layer is closer to the arterial wall. Moreover, Figure 9 (e) 
shows the role of convective flow in drug delivery to (at the reference point) 


the arterial wall tissue. It is clear that convection contributes significantly to 
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«10° 


2 


B 


(a) Drug release, 0 second (b) Drug release, 25 second 


(c) Drug release, 100 second (d) Drug release, 300 second 


Figure 8: Dissemination for different times 
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drug delivery, which must be taken into account in order to achieve better 
results. 

The behavior of mean drug concentration in the inner and outer layers 
of the stent for different diffusion coefficients is shown in Figure 10. This 
figure shows that the reduction of the diffusion coefficient in the inner and 
outer layers of the stent makes the drug enter the arterial wall tissue with a 
delay from the stent. This enables the stent designers to deliver the required 
amount of drug to the arterial wall tissue. 

The average free drug in arterial wall tissue is depicted in Figure 11. 
This figure shows that the drug concentration tends to increase from zero to 
maximum and diminishes afterward. It is interesting to note that in the case 
of fast release of the inner layer (D, = 1 x 10~%em?/s) and slow release of the 
outer layer (Dz = 1 x 10~'cm?/s), this maximum is higher than in the case 
of slow release of the inner layer (D, = 1 x 10~'cm?/s) and slow release of 
the outer layer (Dz = 1 x 107!4cem?/s). 

The arterial wall tissue has a viscoelastic property, which significantly 
affects drug delivery. As a result, we presented this effect in Figure 12. It 
can be seen that the higher the viscoelastic diffusion coefficient, D,, the less 
drug is observed in the arterial wall tissue. In other words, the tissue of the 
arterial wall has a barrier property against drug delivery. 

In order to validate the application of the model presented in this paper, 
in Figure 13, we compare the numerical results obtained for the free drug 
concentration in the arterial wall with the in-vivo experimental data (Sari- 
fuddin and Mandal [24], Tzafriri et al. [26], and Mongrain et al. [18]). It 
is clear that considering the viscoelastic properties of the arterial wall can 
affect the transfer of free drugs from the stent coating layers to the arterial 
wall. These results indicate that the proposed model is congruent with the 
models proposed in Sarifuddin, Tzafriri,Mongrain, although they did not take 
into account the viscoelastic properties of the arterial wall. In other words, 
the model presented in this paper is more consistent with the experimental 
results than other models since the viscoelastic property, which is one of the 
mechanical properties of the arterial wall tissue that acts as a barrier, has 
been considered. Figure 13 illustrates the advantages of this article compared 


to the previous ones. 
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The effect of the diffusion coefficient of the layers of stent coating on free and 


bound drug concentration and the effect of convection on free drug concentration in the 


arterial wall tissue. 
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Figure 10: The effect of diffusion coefficients on the drug in coating layers 
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Figure 11: The average concentration of free drug in arterial wall tissue (the effect of 


diffusion coefficients on drug in coating layers) 
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Figure 12: The influence of D, on the concentration of free drug in arterial wall tissue 
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delivery from stent coating layers into the arterial wall 
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6 Discussion 


In order to discuss the transfer of drugs from a half-embedded two-layer coat- 
ing stent into arterial walls with viscoelastic properties, the following points 
can be specified. In Figure 5, the independence of numerical results from the 
mesh size is discussed. The effect of the velocity vectors on the distribution 
and retention of the drug is shown in Figure 6, which corroborates the results 
in [28]. The effect of filtration velocity on the average drug concentration in 
arterial wall tissue is shown in Figure 7. An increase in time indicates a 
rise in the concentration of the drug in the tissue of the arterial wall, which 
can be seen in Figure 8. As Figure 9 demonstrates, the concentration of the 
bound drug increased and reached a steady state, while the concentration of 
the free drug went up for some time but then decreased. In addition, faster 
release results in higher concentrations of free and bound drugs in the arterial 
wall. Because the inner layer is closer to the arterial wall, the rapid diffusion 
of the inner layer is greater than the concentration of the outer layer in the 
arterial wall tissue. Additionally, convection significantly contributes to drug 
delivery. Figure 10 shows that the reduction of the diffusion coefficient in 
the inner and outer layers of the stent makes the drug enter the arterial wall 
tissue with a delay, which allows the stent designers to determine the amount 
of drugs required for the arterial wall tissue. When the inner layer is faster 
than the outer layer, compared with slow release in both layers, the maxi- 
mum drug concentration is higher, which is observed in Figure 11. Figure 
12 illustrates that the viscoelastic property of arterial wall tissue acts as a 
barrier against drug transmission, which is not discussed in similar works. Fi- 
nally, a comparison between the experimental data and the proposed model 
for drug delivery from the stent coating layers to the arterial wall, which 


exhibits greater agreement than other models, is presented in Figure 13. 


7 Conclusion 


Considering the mechanical properties of the arterial wall due to atheroscle- 


rotic changes provided a deeper understanding of the stent’s performance. It 
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enabled a better agreement between the mathematical model and the experi- 
mental results. A mathematical model is presented to consider the transport 
of drugs from the stent layers to the arterial wall and describes the convective 
flow due to the flow of interstitial fluid in the arterial wall tissue as well as the 
viscoelastic property of the arterial wall. Presenting an upper bound for the 
energy estimate determines the uniqueness of the solution and its stability for 
a limited time. Finally, the numerical simulation verifies the agreement be- 


tween the mathematical model, numerical results, and experimental results. 
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